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We formulate a new analytical method for regularizing the self-force acting on a parti- 
cle of small mass fi orbiting a black hole of mass AI, where fi <^ M. At first order in fj,, 
the geometry is perturbed and the motion of the particle is affected by its self-force. The 
self-force, however, diverges at the location of the particle, and hence should be regularized. 
It is known that a properly regularized self-force is given by the tail part (or the -R-part) 
of the self-field, obtained by subtracting the direct part (or the 5-part) from the full self- 
field. The most successful method of regularization proposed so far relies on the spherical 
harmonic decomposition of the self-force, the so-called mode-sum regularization or mode 
decomposition regularization. However, except for some special orbits, no systematic ana- 
lytical method for computing the regularized self-force has been constructed. In this paper, 
utilizing a new decomposition of the retarded Green function in the frequency domain, we 
formulate a systematic method for the computation of the self-force in the time domain. Our 
method relies on the post-Newtonian (PN) expansion, but the order of the expansion can 
be arbitrarily high. To demonstrate the essence of our method, in this paper, we focus on a 
scalar charged particle on the Schwarzschild background. Generalization to the gravitational 
case is straightforward, except for some subtle issues related with the choice of gauge (which 
exists irrespective of regularization methods). 

§1. Introduction 

We are now at the dawn of gravitational wave astronomy /astropiiysics. Tlie 
interferometric gravitational wave detectors LIGO/^ TAMASOO^^ and GEO600^^ are 
currently in the early stage of their operations, and VIRGO'^^ is expected to be in 
operation soon. Furthermore, R&D studies of a space-based interferometer project, 
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LISA,^) are in rapid progress. It is expected that these interferometers, after their 
final sensitivity goals are achieved, will detect gravitational waves from compact star 
binaries and compact stars orbiting super-massive black holes. 

To fully utilize the information contained in the observed gravitational wave 
data, and particularly for the purpose of precision testing general relativity, it is 
essential to have accurate theoretical predictions of the waveforms.^^ For nearly 
equal mass binaries, the standard PN approximation is a powerful tool to compute 
the waveforms. An alternative method of computing the waveforms is the black hole 
perturbation approach. ^^^^^^ This approach is very effective, in particular, when the 
mass ratio of the objects composing a binary is extreme. In this paper, considering 
the case of such extreme mass ratios, we propose a new method for calculating the 
corrections to the force acting on the small mass body (which is treated as a point 
particle) that are induced by the field generated by the particle itself, the so-called 
self-force corrections. 

In the black hole perturbation approach, one can appeal to the energy- angular 
momentum balance argument to evaluate the radiation reaction to the orbit of a 
particle, namely, by equating the rates of change of the energy and angular momen- 
tum of the particle with those carried away by the gravitational waves emitted by it. 
However, the balance argument is not sufficient. First, the radiation reaction to the 
Carter constant can not be calculated with this method. The Carter constant is the 
third constant of motion of a test particle in the Kerr spacetime. The other two con- 
stants of motion, the energy and angular momentum with respect to the symmetry 
axis, are associated with the Killing vectors of the background spacetime, and their 
rates of change may be evaluated from the gravitational waves emitted to infinity or 
absorbed at the black hole horizon. In contrast, there is no correspondence to any 
such Killing vector in the case of the Carter constant. Hence, its rate of change is 
not directly related to the waves emitted. Second, and most importantly, the bal- 
ance argument can yield only time-averaged rates of change of the two constants of 
motion, while there are many situations in which knowledge of the actual radiation 
reaction force, as well as the so-called conservative part of the self-force, at each 
instant of the orbital motion becomes necessary. 

Consider a particle having either a scalar, electromagnetic or gravitational charge. 
The orbital motion of the particle will create a field at first order in its charge, and 
the motion will be affected by its self-field. This self-field, however, is divergent 
at the location of the particle. Hence, the force due to this self-field is apparently 
ill-defined. It is known that the self-force in the vicinity of the particle may be 
decomposed into the so-called direct part and the tail part, and that the correctly 
regularized self-force is given by the tail part. The justification of this prescription 
is given in Ref. 14) for the scalar and electromagnetic cases, and in Ref. 15) and 16) 
for the gravitational case. 

In the scalar case, the Klein-Gordon equation is hyperbolic from the very begin- 
ning. In the electromagnetic and gravitational cases, the field equations can be put 
into hyperbolic form by choosing the Lorenz gauge (often called the 'harmonic gauge' 
in the gravitational case). In general, for a hyperbolic equation, the retarded Green 
function G^^^{x, x') can be split into two distinct pieces (at least locally when the two 



Self-Force Regularization 



3 



points, namely, the field point x and the source point x' are sufficiently close), which 
is called the Hadamard form. One piece has support only on the future light-cone 
of x', and the other piece has support in the interior of the future light-cone of x' . 
The former gives the direct part and the latter the tail part. Recently, an equivalent 
but more elegant decomposition of the Green function was proposed,^'') in which the 
direct part is replaced by the S'-part and the tail part by the i?-part. The S-part is 
defined by adding a piece that has support outside the light-cone in such a way that 
it does not contribute to the self-force when it is subtracted from the full field. The 
remaining part is called the i?-part, which now has support outside the light-cone as 
well. The advantage of this new decomposition is that the S'-part is symmetric with 
respect to x and x' , and it satisfies the same equation as the retarded Green function. 
This implies that the i?-part now satisfies a source-free, homogeneous equation. 

Thus our task is to evaluate the tail or R-part of the field and the self-force due 
to this. However, we do not have any systematic method to compute the i?-part 
directly. In contrast, there exist several schemes to compute the full retarded field 
in Schwarzschild and Kerr spacetimes.®-*"^'^-*'^^-*"^^-* Therefore, what one can do is 
subtract the direct or 5-part from the full field to compute the regularized self-force. 
There have been many investigation of this method. 

Because both the full field and the S'-part diverge at the location of the particle, 
it is necessary to develop a regularization scheme to compute the difference between 
the two. The most successful scheme of regularization that has emerged is the mode 
decomposition (or mode-sum) reg 

ularization.22)'23) The full field can be decomposed 
into partial waves by using the spherical harmonics y^„j(J7). The contribution to the 
force from each i-mode does not diverge in the coincidence limit of the field point 
with the location of the particle. Now, if we subtract the S-part from the full field, 
before summing over i, the divergence disappears. Hence, we can perform summation 
over i to obtain the expression for the regularized force. 

The S-part can be calculated only in the form of a local expansion, that is, at 
field points sufficiently close to a point on the orbit, and therefore it is necessary to 
extend it over a sphere containing the orbital point to obtain the spherical harmonic 
coefficients. In recent years, a method to carry out its harmonic decomposition has 
been developed in the Schwarzschild case,^^-'"^''-* and it was recently extended to the 
case of the Kerr background.^^) The full field is calculated using either the Regge- 
Wheeler-Zcrilli or Teukolsky formalism, which uses the spherical (or spheroidal in 
the Kerr case) harmonic decomposition. 

To this time, the regularized self-force has only been calculated numerically. 
These numerical results are important. However, an analytical understanding will 
be very useful. Analytical results obtained to this time are restricted to the case of 
particles orbiting in very special orbits (mainly circular or radial infall), and there 
axe no results for the general case. This is primarily due to the mismatch in the 
schemes used for evaluating the full Green function and the S'-part. The Regge- 
Wheclcr-Zcrilli and Teukolsky formalism rely heavily on the Fourier decomposition 
of the time-dependence by taking full advantage of the stationarity of the back- 
ground spacetime, and therefore the full field is calculated in the frequency domain. 
Contrastingly, the S'-part is in the time domain. If the orbit is specified a priori, it 
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is in principle possible to obtain the full field in the time domain by explicitly per- 
forming the integration (or summation) over the frequency oj. However, in practice, 
the explicit integration over w is possible only in very special cases, such as circular 
orbits.29).30) 

In this paper, we propose a new method to go from the frequency domain to 
time domain, and regularize analytically in order to calculate the self-force com- 
pletely analytically. Though we employ the PN expansion method, it is formulated 
systematically in such a way that the order of the expansion can be taken arbitrarily 
high, as long as it is kept finite. In fact, our method is more effective in the far zone, 
and thus it supplements the previous numerical method, which is effective near the 
plunge. For simplicity, we restrict our analysis to the case of a geodesic orbit in the 
background spacetime, but removing this restriction is straightforward (although, 
of course, the equations become much more lengthy). Furthermore, we focus on 
the case of a scalar charge in order to avoid the gauge problem. There is a subtle 
problem associated with the choice of gauge in the gravitational case, 3^^' ^^-^ but our 
method for integration over oo is equally applicable to the gravitational case, despite 
this problem. 

The key idea of the new approach developed in this paper is to separate the 
retarded Green function in the frequency domain into two distinct pieces, in analogy 
to the S-R decomposition in the space-time domain. We call them the S'-part and 
the -R-part. The former contains all the singular terms to be subtracted, while the 
latter satisfies the source-free, homogeneous equation. In particular, once the PN 
order is specified, the contributions from only a finite number of £ is necessary to 
evaluate the -R-part. The most important point of this new decomposition is that 
the S'-part in the frequency domain is given in the form of a simple Taylor series with 
respect to lv multiplied by exp[— zu;(t — t')]. Therefore, the integration over oj can 
be performed easily for such terms. They just produce 5{t — t') and its derivatives. 
Using this technique, we can obtain the 5-part in the time domain relatively easily. 
Then, the regularization is done by subtracting the 5-part from the thus obtained 
S-part in the time domain. 

This paper is organized as follows. In 32 we describe our new decomposition of 
the retarded Green function in the frequency domain. In ^ focusing on the scalar 
case, we demonstrate our regularization method. We first integrate the force due to 
the S'-part over lo to obtain the harmonic modes of the force in the time domain, 
then subtract the S-part mode-by-mode, and finally sum over d. to obtain the force 
due to the (5 — S)-part. We do not discuss the i?-part, because it is finite from 
the beginning. The final section, ^ is devoted to conclusions and discussion. Some 
formulas and proofs of several propositions used in the text are given in appendices. 
For readers' convenience, formulas for 4PN order calculations can be found at the 
ht tp : / / www2 . yukawa. kyoto-u . ac.jp/~misao/BHPC/. 
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§2. New decomposition of the Green function in the frequency domain 

We consider a point scalar charge q moving in the Schwarzschild background 



2M\ 



r I 



l- — \di^+[l- — ] dr^ + r"^ {d9^ + sin^ edif'^) , (2-1 



r 7 



where {x"} = {t, r, 9, cp} are the Schwarzschild coordinates, and M is the black 
hole mass. The full scalar field induced by this charged particle is given, using the 
retarded Green function, by 



= -q / dTG*^"(x,z(r)), 



(2-2) 



where r is the proper time of the particle, and G^"'^(x, x') satisfies the Klein-Gordon 
equation. 



5(4) (x-x') 



-9 



(2-3) 



with retarded boundary conditions. The full Green function is represented in terms 
of the Fourier-harmonic decomposition as 

Gf"^x,x') = J 5] 5rjr,r')y,„.(0, 0)^4(0', <^') . (2-4) 

£m 

Here, Yim{0,(l)) are the ordinary spherical harmonics. Then, Eq. (|2-3j) reduces to an 
ordinary differential equation for the radial Green function as 



/ 2M\ d^ 2(r -M) d 
\ r J dr"^ dr 



( 



UJ 



£{£+!) 



V 



2M 



--<5(r-r'). (2-5) 



The radial part of the full Green function can be expressed in terms of homo- 
geneous solutions of Ea. ()2-5() . which can be obtained using a systematic analytic 
method developed in Ref. 18). We have 



-1 



W (^-^ ^\ (C(r)<p(r')^(r' -r) + ct>''n^{rWn{r')e{r - r')) , 



d 



(2-6) 



Here, the in-going and up-going homogeneous solutions are denoted, respectively, by 
(f))^^ and (/>^p, and v is called the 'renormalized angular momentum', ^^)' '^^^ which is 
equal to i in the limit Mto — > 0. 
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We express the homogeneous solutions (p^^ and (p'^^ in terms of the Coulomb 
wave functions cj)'^ and </>^^~^ .13), 18) 

<p = 7.0c+'J.0r"'- (2-7) 

The properties and the relations of the coefficients {a^, (3^, Iv, 5y} are studied ex- 
tensively in Ref. 18) and 19). (The function (j)'^ is denoted by i?^ in Ref. 18)) Here, 
we need to stress a remarkable property of the wave function which becomes 
manifest when we consider the PN expansion, i.e., when cf)^ is expanded in terms 
of z := ojr and e := 2Moj, assuming they are small [0{v) and 0(t>^), respectively]. 
In the expression of its PN expansion, := {2z)~'^(j)'^ contains only terms that are 
integer powers of z and e, and there are no terms like log z. In fact, we find that 
this condition, that logz is absent, uniquely specifies a single solution of the radial 
homogeneous equation. As is explained in Appendix this fact can be used to 
compute simultaneously determining the eigenvalue i^. Furthermore, this PN 
expansion turns out to be a double Taylor series expansion in and e/z, i.e., with 
only positive powers of Lo^. With the normalization ^ 1, at the leading order, 
the expansion is 



2 ,4 



z^ -M- 10)ez 



^ ^ 2(2£ + 3) 2z ^ 8(2£ + 3)(2£ + 5) ^ 4(2£ + 3)(£ + l) ^ 



15£2 + 15£-11 2 



2(2£- 1)(2^ + l)(2£ + 3) 



e^ + ---. (2-8) 



The solution (j)^'^ ^ can be obtained through the replacement £ —£ — 1. The 
Wronskian of (j)'^ and (p^'^^^ becomes 

Of _|_ 1 

(496^6 + 1488^5 + 1336^^ + 192^^ _ 757^2 _ q^^^ ^ 333)^2 
^ 16(2£-l)2(2£ + l)(2£ + 3)2 ^ ' 

We note, however, that the general expression for the PN expansion of Pc'^^^ with 
this method (i.e., requiring the absence of log z terms) becomes invalid at the (£— l)th 
PN order. For small i (< PN + 1), the computation must be done following the 
systematic method given in Ref. 18). The respective results for the i = and 1 cases 
are 

^„ 7 7^2 7e 7z'-^ Uez 

= \ \ 

9 54 27z 1080 27 
^ = ]o(z) + IH \ \ , 

7 2 .r. 49 23263 o 

u = — e^ + ---, ioW = + e^ + ---, 2-10 

6 ' 162 29169 ' ^ ^ 

and 
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-Z/-1 



30^3 



29^2 



19 o 3 117443 



z^ = l-7:^e^H , = -- + ———e^ + 



30' ' ' -- 2 ' 216600 

Here, jo (-?) and ji (z) are the spherical Bessel functions. 
We now divide the Green function into two parts, as 



1^ 



(2-11) 



(2-12) 



where 



-1 



€{r)rc{r')e{r' - r) + 4>-''-\r)(t>l{r')e{r - r' 



(2-13) 
(2-14) 



Here, we have assumed that ^ and ^ 0, and we have also introduced the 

coefficients {j3i,-, ^u} '■= {Pu/ctu, Ivl^v)- Using results obtained in Ref. 18), we obtain 
the behavior of the coefficients in the PN expansion as 



The functions 0c ^^^^ respectively, of 0(f ) and 0(u 



(2-15) 

(except at 

£ = 0). Therefore, the three terms in the i?-part of the Green function become, 
respectively, of 0(t;®^), 0{v'^^^^) and 0(f^^"'"^) relative to the S-part. 

The part that we need to consider in the regularization of the force is just S. 
There is no divergence associated with the remaining i?-part, and it terminates at 
finite i as long as we restrict our consideration to finite PN order. Moreover, the 
^-part satisfies the homogeneous radial equation. This fact will be an additional 
advantage of the present method when we consider extension to the case of grav- 
ity. Because the i?-part is a homogeneous solution, we can apply Chrzanowski's 
method^''-* to reconstruct the metric perturbations even in the frequency domain. 
We discuss this point in more detail in a separate paper. 

§3. Computation of the S-part 

We now compute the force due to the 5-part for a general orbit. When expanded 
in terms of the spherical harmonics, the force corresponding to the iS-part (the S- 
force) is known to take the form^^^ 

lim F^i, 



x—*zo 



(3-1) 



where is the ^-mode of the 5'-force, L = £+1/2, and Aa and 3^ are independent 
of L. When summed over £, the A-term gives rise to a quadratic divergence, and the 
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B-ieiTo. diverges linearly. For large I, Dae is at most of 0{L~'^), and hence it contains 
no divergence. Because the 5-part can be calculated only locally, its extension to the 
entire sphere involves some ambiguity. As a result, the coefficient of each i-mode, 
Dai, depends on the method of extension, but the final result after summation over 
i, which is determined by the local behavior of the field near the source location, 
does not. It is known that 

oo 

J2Da,i = 0. (3-2) 

e=o 

The difference between the S'-force and the S-iorce should be finite, because the R- 
force is finite. Thus, in general, the S'-force must take the same form as the S'-force 

lim F^^ = AaL + Ba + Da,i ■ (3-3) 

Below we confirm explicitly that both Aq, and for the S-force coincide with those 
for the S-force. Therefore, the force due to the S-part minus the S-part, which is 
finite, is given by 

TPS-S 



1=0 e=o 

3.1. Force in the time domain 

To obtain an expression for the S-force in the form of Eq. (|3-3|) , it is necessary 
to perform the u> integration explicitly. Here, the key fact is that there appears no 
fractional power of uj in the S-part. This is because we have chosen (p'^ and 4'^'^^^ as 
the two independent basis functions. As noted above, except for the overall fractional 
powers z'^ and z~'^~^, they contain only the terms with positive integer powers of 
(j^. When we consider a product of these two functions, uj contained in the overall 
factors z'^ and z~'^~^ just produces uj~^, which is canceled by uj from the inverse of 
the Wronskian. Thus, g£rnuj{r,r') is expanded as 

oo 

9l^{r,r') = Y,^^'Gimk{r,r'). (3-5) 

fc=0 

Therefore, the integration over uj can be performed easily, as we shall show now 
explicitly. The Fourier transform of simply produces 

dww^^e-^'^^*-*') = 27r(-l)"a|"<5(t - t'). 

Differentiation of the delta function in the expression above can be integrated by 
parts to act on the source term. Thus, we can express the S-force in the time 
domain as 



m,k 

(3-6) 
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Here we have inserted the projection tensor P^^ = Sa^ + UaU^, where is the 
four- velocity, so that the normahzation u^Ua = — 1 is maintained. Now that we have 
the 5- force given in the time domain, the meaning of the coincidence hmit r — > z'' (t) 
is transparent. 

Here one can set z^{t) = 7r/2 without loss of generality. Then, nothing that 
d^Y(m{iT/2,(p) = imYem{'^/'^,f) and 9^Y^^(7r/2, = -imY^^(7r/2, the summa- 
tion over m can be done by using the formulas 

y Y,^ ^)" = l { ="7f ' 

^ £m ^ 2 ' ^7 I 0, for J =odd, 

m=—£ 

e 

y deY,^ {9, (^)[^^^^ Y,*^ g, ^) = 0, (3-7) 

m=—£ 

where A(„) {£) is a polynomial function of i of order n + 1 defined by 

n=0 ^ 

3.2. Separation of the A-term 

Before performing the operation discussed in the preceding subsection, we sepa- 
rate the A-term from the other contributions. This can be done easily by using the 
fact that only the A-term has a jump in its values at the coincidence limit, depending 
on the direction from which the source point is approached. We divide the S-part 
Green function into symmetric and anti-symmetric parts as 

gLjr, r') = g^^Jir, r') + sgn(r - r')g^^J{r, /), (3-9) 
where sgn(y) = ±1 for y ^ 0, and 

aUtHry) = u!^-u-u mr)<f>-^-Hr') ± <Py-\rWAr')] • (3-10) 

Then the force due to the anti-symmetric part must coincide exactly with the A- 
term. We know that the A-term for the S'-part has a simple form proportional to 
L, while the expression for the S- force given in terms of the Green function is more 
complicated. The reason why such a simple result for the A-term is recovered is 
explained in Appendix ^ We concentrate now on the symmetric part, which is 
responsible for the B and D-terms. 

3.3. The {S — S)-part of the force 

The result for the S'(-|-)-part of the force is 

2 J. oo 2 r r 

n=0 n=0 

and 

'^'^ " n-(l-2M/ro) u^rf^^' ' ^ ' 
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where the coefficients K^^^, whose upper index (i) represents the PN order, are 
formally given by 

E di^^^\l) {5er UK (3-13) 

1 M 

Here, the quantities da are some functions of i, vq = z^Uq), = 1 — —7j, U = — , and 

ro 

£ and C are, respectively, the energy and the angular momentum of the particle. To 
obtain these expressions, we have used the first integrals of the geodesic equations, 

dz''it)Y_f^ 2M y 1 / 2My/ _^ £2 



dt J \ z'-{t) J E"^ \ z'-it) J V ^''(t) 



fz'-(t)2V -2''(i)y' dr l-2M/z^{t)'' 
and we have also reduced higher-order derivatives with respect to t by using the 
equations of motion 

d'^z'it) _ 2M '^^ \ + (i 



dt^ z^'itf V z'-it) J £'^z'-{tf \ z^it) 
3M / 2M \ V 



As long as we ignore corrections that are higher order in /i, we can assume that the 
orbit is momentarily geodesic. Here, we find that the 5-part force is written solely 
in terms of the orbit at the location of the particle (that is, no tails!). Hence, for 
the correction at lowest order in /i, the force coming from the {S — 5)-part is also 
written in terms of the position and the velocity of the particle. We do not have any 
terms with a positive power of i, as expected. However, this cancellation looks rather 
miraculous in the present formulation. In AppendixEl we give a brief explanation of 
why the terms with positive powers of i are absent. From the asymptotic behavior 
for large i, we can read off the coefficients Ba- These coefficients are identical to 
the results obtained previously. 22)'25) As mentioned earlier, a separate treatment is 
necessary for small i. To summarize, the {S — 5)-part of the force is given by 

CO CO / / 1 \ \ 

^-"'' = Ei™ {F!/-Fle)=Y.[^le-^aU+^)-Bj 

CO 

= E(<i"^-^") (3-16) 

e=o 

The summation over i is performed by using the decomposition into partial 
fractions and the formulas 

°° 1 R °° 1 R 

^_^^ = ^ " M!' E (7317^ = 2 (2 -1)^ (3-17) 
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where n is a positive integer, and is the Bernoulh number defined by 

+ --l = y 22"(22- - 1)7^x2". (3-18) 
e^-1 2 ^ ^ \2n)\ ^ ' 

n=l ^ ' 

The reason why odd powers of I do not arise is as follows. If we have a term hke 

OO ^ 

then, due to the symmetry under £ ^ — + 1), we also have a term 

OO ^ 

Jj-£-l + fe/2)"' 

where /c is an integer. If n is odd, these two contributions combine to leave a 
summation of finite terms. Hence, there remains no infinite summation of odd power 
terms in the final expression. 

The result for the {S — S')-part of the scalar self-force is 

" n=0 " n=0 



and 

Fr~^ = TT-Avy^Ff''^ - -^^<f (3-20) 



where the coefficients '^^'^ whose upper index, i, represents the PN order, are 
formally given by 

j+j+fc=n \ / 

CS-S(n)^ ^ efe-^)(5£r f^^jV^ (3-21) 



Here are some constant numbers. 

Once we obtain the general expression for the {S — S')-part of the force, compu- 
tation of the remaining _R-part is rather easy, because only terms up to a finite value 
of i contribute to the force for a given PN order. Then, the i?-part of the force, 
which is what we want in the end, is given by 

= Ft' + F^. (3-22) 



Now we discuss a technical but important property of the iS-part of the Green 
function. The two independent radial functions at leading order in e are given by the 
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spherical Bessel functions ji{z) and ni{z). Both (j)'^ and (t>c'^^'^ are given by linear 
combinations of these two independent solutions. Up to 0{e^), we have (f)^ oc ji{z), 
while 4>c'^^'^ oc e^'^je,{z) + Cini{z), where is a constant of order unity. Here, in 
passing, we note that the leading term of (t)c^~^ the PN expansion comes not 
from the term e~^je,{z) but from ni{z) for I > 2. This is because the ratio of the 
two terms is ji{z)/{en£{z)) oc = 0{v'^^~'^) in the PN expansion. At first glance, 

the existence of the term e~^je{z) seems problematic, since it would naively lead to 
a term in the Green function that behaves as 1/M. Collecting the leading terms in 
e, we find that the Green function has a term proportional to 1/M of the form 

This expression is identical the radiative Green function in Minkowski space, except 
for an additional multiplicative factor 1/{Muj). Therefore, after summation over I 
and m and integration over oj, we find that this part is 



oc 



dt' 



(5(t - i' - |x - x'l) 5{t-t' + \yL 



|x — x'l |x — x'l 

e{t - f + |x - x'l) - e{t - f - |x - x'l 

|x — x'l 



At lowest order in M, the trajectory of a particle is a straight line in the Minkowski 
background. Because the above expression for the leading part of the Green function 
is Lorentz invariant, we can choose this straight line to be static, without loss of 
generality. Then, it is easy to see that the component of the field proportional to 
1/M is constant. Hence, this part does not contribute to the force. 

Another important property of the force due to the {S — S')-part is that it 
contains only the conservative part of the force. To show this, we use the fact 
that the equations of motion take the form if = (a function even in vl') and £ = 
£ = 0, at leading order in /i [see Eqs. (|3-14p and ()3-15() ]. Here, the dot represents 
differentiation with respect to t. Recalling the general expression for the S-force in 
Eq. (|3-6p . we see that the t-component contains an odd number of time derivatives, 
Q2k+i I Q^2k+i ^ the r-component an even number of time derivatives plus one radial 
derivative, d"^^^^ / {dt"^^ dr) , and the (/^-component an even number of time derivatives 
plus one (f derivative, 8"^^+^ /{dt^^dif). Now, using the equations of motion, the t 
derivative can be replaced by 

d d d d 

!L = z'-{t)^ + -f{t)^ + z^{t) 



z^it)^+z^{t)^-zmi^t), 



and the if derivative by 



d 

— = -\-im . 
dip 
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Notice that one differentiation with respect to i or 93 changes the total power of 
and m by an odd number, while a differentiation with respect to r does not. Notice 
also that only the terms even in m remain after summation over m. Therefore the 
£-mode of the S'-force takes the form 

Ff, = Tti{r,£,C)u'', Ff, = J^re{r,£,C), F^, = J'^e{r,£,C)u' . (3-23) 

The S'-part of the force is known to have exactly the same form. This implies that 
the {S — 5)-part of the force also takes the same form. Thus, after summing over i, 
we conclude that the final form of the {S — S')-part of the force is 

F/-^ = .F^(r,f,£)^/^ F,^-^ = .F,(r, £), = .F^(r, £) (3-24) 

We can now explicitly show that the above form of the force implies the ab- 
sence of a dissipative reaction effect. In other words, the force is conservative. The 
equations of motion to 0(/u^) are given by 

fj^^ut' = F^, (3-25) 
ar 

where is the perturbed four velocity and D/dr is the covariant derivative. Then, 
we obtain the evolution equation for the perturbed energy £ := —^t^u^ as 

f = -f.^O") = = -^*)|. (3-26) 

where = {dtj^ is the time-like Killing vector. This equation is integrated to give 

£ = £- j Tt{r)dr. (3-27) 

Here, £ is an integration constant, which we can interpret as the unperturbed energy. 
In the same manner, for the perturbed angular momentum £, we obtain 

/r 
J^^{r)dr. (3-28) 

Thus we find that there is no cumulative effect on the evolution of the energy and an- 
gular momentum of the particle. In other words, a force of the form (|3-24|) preserves 
the presence of the constants of motion £ and C. Concerning the radial motion, 
u'^ can be expressed in terms of fiut = £ and /i = £ by using the normalization 
condition of the four velocity. We have 



~r 
/U U 



.2 .0. ox 1/2 



: £^-{l~ 2M/r) ( 1 + C'/r' ) . (3-29) 
Thus, vl' is obtained function of r. 
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§4. Conclusions and discussion 

The present work is an attempt to make progress toward a more realistic (calcu- 
lation effective) analytic scheme for constructing orbits, taking into account radiation 
reaction effects. The key idea proposed in this paper is a new decomposition of the 
Green function into S and i?-parts, given in Eq. (|2T2)) . This new decomposition 
relies on a systematic analytic approach to the black hole perturbation developed in 
Refs. 13), 18) and 19). The new decomposition is not identical to the usual S and 
R decomposition,^'^) but they share certain properties. The S-part is singular and 
symmetric, and it satisfies the same inhomogeneous equation as the Green function. 
The i?-part is regular, and it satisfies the source-free equation. Considering a scalar 
charged particle, we showed that the S-part of the self-force can be evaluated analyt- 
ically in the time domain and that it yields the same regularization parameters 
and Ba in the mode-decomposition regularization as the usual 5-part. This implies 
that the S-part contains all the singular behavior of the original S'-part. Also, we 
showed that the self-force due to the {S — 5)-part is conservative. Moreover, we 
found that the i?-part of the force valid up to the {(. + 0.5)th Post-Newtonian (PN) 
order can be obtained by taking account of only the spherical harmonic modes up 
to £-th order. 

The analysis here is restricted to the self-force due to a scalar field for its sim- 
plicity. The extension to the electromagnetic case is straightforward, although the 
computation of the force from the master variable of the perturbations^'^) becomes 
more tedious. The gravitational case, however, is more complicated, because in that 
case, the self-force is gauge dependent. In the scalar case, as is manifest from our 
calculation, the computation of the S'-part is in fact unnecessary if we make use of 
the fact that the non-singular part of the regularization parameter Dai for the S'-part 
vanishes after summing over I modes. By virtue of this fact, the simple prescription 
of subtracting A- and S-terms is valid in this case. In the case of gravity, before sub- 
tracting the 5-part from the S'-part of the force, we need to adjust the gauges that 
are originally different. A natural prescription for this is to transform the S'-part, 
originally given in Lorenz gauge, to the gauge in which the S-part is computed. Due 
to this gauge transformation, the parameter D^i as well as the other regularization 
parameters will be altered. Then the contribution to the self-force from Dan will not 
vanish in general. This needs to be studied in more detail. 

Related to the gauge dependence of the gravitational self-force, there arises a 
conceptual question: What kind of gauge-independent concepts are contained in 
this otherwise gauge-dependent quantity? There are several constants of motion for 
geodesies in the background black hole spacetime. These constants of motion evolve 
after we incorporate the self-force. However, the secular change of "the constants 
of motion" has a gauge-invariant meaning. It can be evaluated by merely taking 
account of the force due to the i?-part, or we can simply use the radiative Green 
function. In addition to the information regarding the constants of motion, the 
self-force may contain other gauge-invariant information. This question needs to be 
answered, although it is not a problem specific to our present method. We may find 
that the physical information contained in the gauge- dependent self-force by itself is 
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very limited. Even if this is the case, calculation of the self-force will be a necessary 
step to develop a black hole perturbation theory to second order in the mass of the 
orbiting particle. 

One of the main advantages of our method is that it allows a successful im- 
plementation of a systematic post-Newtonian expansion technique in the black hole 
perturbation. There could be criticism of our approach in regard to the limita- 
tion of the PN expansion itself. The black hole perturbation is considered to be 
a method complimentary to the standard post-Newtonian approximation. In this 
sense, one might think that there is no point in using the PN expansion in the black 
hole perturbation approach. However, we should stress that the PN expansion of 
the perturbation in the black hole spacetime can be systematically extended to an 
arbitrarily high order without any conceptual difficulties. Hence, if we can make 
use of this advantage, problems far beyond the validity of the finite order standard 
post-Newtonian approximation can be investigated. 

In an actual computation, the achievable PN order will be limited. Then, the 
question will be the speed of the convergence of the PN expansion. In this regard, 
there are a couple of encouraging pieces of evidence. First, we mention the issue 
of the innermost stable circular orbit (ISCO). At present, the calculation of the 
orbital frequency at ISCO up to 3PN order has been madc.^^^ The result is in rather 
good agreement with the numerical result for an equal-mass binary. Second, we 
mention the energy loss rate of a particle in a circular orbit in the Schwarzschild 
black hole evaluated from the asymptotic waveform. Although the convergence 
of the PN expansion becomes slower and slower for a smaller orbital radius, there 
is no evidence of the failure converge, even at ISC0.36) Hence, if we can develop a 
systematic method of evaluating the gravitational self-force in the PN expansion, its 
range of validity should be very wide. 
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Appendix A 

An Easy Way to Find (jf for Sufficiently Large I 



Consider a radial function cjf {z) which satisfies 



2M 



^ 2(r_ 



M) d 



dr 



+ 



1 - 



2M 



) 



6*^(0 = 0. (A-1) 



As long as we consider sufficiently large this radial wave function can be completely 
specified up to an overall normalization by the requirement that (p^^ := {2z)^^(f)'^^ 
docs not contain logz in its PN power series expansion with respect to and e/z. 
The same condition simultaneously determines the renormalized angular momentum 
u i. In fact, the equation for becomes 



[z^dl + 2z{v + l)d,] 



l(2-l).^a.^ + i{(4. + 3)-i(2. + l)} 



zdz 



+{e-iy){e + u + i)(i--)-z^ + -(i--) 



(A-2) 



The right-hand side of this equation is of higher order in the PN expansion. Substi- 
tuting the Taylor expansion oi^^ with respect to e/z and z^ into the above equation, 
the coefficients are determined order by order. However, the z-independent terms 
in vanish on the left-hand side. Therefore, the terms which become zeroth order 
in z should also vanish on the right-hand side. This condition determines v order 
by order. If, however, i is not sufficiently large, terms proportional to z"^^"^ arise. 
For such terms, the left-hand side is suppressed by a factor of O(e^). Hence, the 
iteration scheme becomes invalid for small and we need to go back to the original 
method presented in Ref. 18). 

Appendix B 

A-Term 



Here, we consider the ^-term extracted from the 5-part. As we have explained 
in the main text, the A-tcrm corresponds to a jump of the field. Therefore, it is 
given by the antisymmetrized Green function 



2WemUrcAc) 

We introduce a new function, Xc '■= f(l>c, which satisfies 

[52.+a;2-y(r)] Xc(r*) = 0, 



^ mr)cp-'-Hr') - c^-'-'irmr')] . (B-l) 



-H) 2(r - M) 



(B-2) 
(B.3) 
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where r* = r + 2M ln(r/2M — 1). The Wronskian is written in terms of r* and Xc as 



WemU€,'Pc''-') = ( ^Xc''-\r*))x'An - i^Xcin ]X-X\n- (B-4) 



We expand the function g^^J ij^-, in a power series with respect to r* — r*' as 

\r,r')=Y,9n{r'W-r*')\ (B-5) 



Imu) 

n>0 



5„(r')=la,",5fi:^(r,r') 



(B-6) 



The higher-order derivatives with respect to 9^* in Qui'''') can be reduced by using 
Eq. ()B-2|) . Hence, either one or zero r*-derivatives remain in the end. Setting r = r', 
the terms with no r* derivative vanish, while the terms with a single derivative yield 
the Wronskian, which cancels the denominator. As a result, we obtain a rather 
simple expression for gf„. In fact, we have 

1 

9i = 



2r'2 ' 

2M\ 1 

92 = 1 



1 



fl'3 



12r 



i-uj^ + V{r'))-2il 



/2 



^'3 



Note that only even positive integer powers of uj appear. Therefore, let us 
consider terms proportional to uj"^^ for a given A^. Because the factor uj"^^ arises 
only from the elimination of the 2A^ derivatives, and because only a single derivative 
can remain at the end of the calculation, u;^^ can be contained only in 9n(^') with 
n > 2N + 1. Conversely, this means that we have 

AT-l 



/ /\ r * *r\2N 2n f l\r * *l\2N 

g2N[r){r -r ) = 2^ w a„(r )(r -r ) , 

52JV+i(rO (r* - r*')'^-*-^ = ^n(r')(r* - r*')'^+^ , (B-7) 

re=0 

where a„(r') and 6„(r') are independent of uj. We replace lo with a time derivative, 
which acts on r' = z'^{t). The force is calculated by differentiating the potential 
once. Hence, the term proportional to g2N vanishes in the coincidence limit, because 
it contains 2A^ — 1 derivatives at most. As for the the term gi27V+i) which contains 
2A^+ 1 derivatives, each derivative must act on each factor of r* — r*' in (7-* _7-*')27V+i 
to give a finite result. Therefore, in gf2Af+i('^')> only need to keep in the 
coincidence limit. Thus, we can simplify the coefficients Qn as 
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Substituting this expression into Eq. (|3-6|) . and summing over the m- modes, the 
anti-symmetric part of the force can be calculated as 



a(. 



n=0 



2M Y / dz''{t)/dt 
\l - 2M/z'{t) 



2n 



ATrz^{t){z-{t) - 2M)(1 + £2/z-(t)2) ' 
where we have used the equation of motion, Eq. (|3-14|) . 

Appendix C 

Absence of Large Powers of d 



(B-9) 



We have z = or and e = 2Muj, which contain iv implicitly. This uj is replaced 
by a time differentiation, which produces mQ^p. That is, lo is effectively of 0{i). In 
other words, we should regard z to be 0{t) and e to be 0{t}. Then, considering the 
radial function given in Eq. (|2-8() . one might think that the final expression for the 
force would have terms with large positive powers of L We explain here why this is 
actually NOT the case, by analyzing the leading I behavior of and 

Let us introduce "^f^ by 



exp 



exp 



druj^" 



(C-1) 



Then the equation for 'F'^ gives 

'£(£ + 1) z2 



+ 



2z-e 



1/2 



z{z — e) {z — e)2 z{z — e) z{z — e) 
Applying the large £ asymptotic expansion to this expression, we have 



(C-2) 



+ 



z{z — e) {z — e)^ 



1/2 



(C-3) 



where % = 0{£°) and Sl^" = {!/£). 

We find that, at the leading order in £, the condition that the term of 0{l/z) in 
should vanish determines in the sense of the PN expansion: 



u = uo + 0{£^), uo = £ 



lei 



+ 



(C-4) 



Also, is given by 



z z 



+ 0{z^ l£^^ 



2£ 8£3 16^5 



+ 
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Substituting this expression into Eq. (|(M|) . the result is 

, ze 



which coincides with Eq. H2-8() in the large (. limit. Another independent solution 
can be obtained through the replacement v — > —u—1, i — > +[• • • ]^/^ — > 

-[•••]l/2 in Eq. (in^ . 

The product of <P^{r) and 'P^'^^^{r'), which appears in the Green function, be- 
comes 

^''{z)<p-''-\z') = expU' dz^^iz) + J' dz' ^o^-^z'U X 0(£°) . (C-7) 



As clearly seen from the expanded form of iP'q given in Eq. ()(>5|) . we have 'Fq\ 
(i/uj) X]^o('^/^)^"^"('^)' where the quantities C„ are functions of r, and \I^q'^~^{z) = 
—^q{z) + 0{\/tj. This implies that, when z and z' are sufficiently close, we have 

f dzW^{z) + r (i/a^o-"-i(z') = £(r - r') Y,{uj/irFn{r,r') + 0{f), (C-8) 

•' n=0 

where the are functions of r and r' that are independent of a; and £ and regular 
in the coincidence limit, r — > r' . Therefore, we obtain 

oo oo 

<P''{z)<p-^-\z') = 5^5^r(r - r'y{uj/irFs,n{ry) + 0{£'), (C-9) 

s=l n=0 

where the Fg^n are functions of r and r' . The terms apparently of 0{i^) are always 
associated with the factor (r — r')'*. Therefore, the time differentiation must act 
on r' = z'^{t') at least s times. Otherwise such terms vanish in the coincidence 
limit. Each time differentiation, however, reduces the power of u by 1, and hence 
it produces the factor l/l. Thus, the terms that appear to be 0{l^) turn out to be 
0{f) in the end. 



1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 

13 
14 



References 

LIGO web page: http://www.ligo.caltech.edu/ 
TAMA300 web page: http://tainago.mtk.nao.ac.jp/ 
GEO600 web page: http://www.geo600.uni-hannover.de/ 
VIRGO web page: http://www.virgo.infn.it/ 
LISA web page: http://lisa.jpl.nasa.gov/ 

C. Cutler et al, Phys. Rev. Lett. 70 (1993), 2984; astro-ph/92 08005| 

L. Blanchet, Living Rev. ReL 5 (2002), 3; gr-qc/0202016 

T. Regge and J. A. Wheeler, Phys. Rev. 108 (1957), 1063. 

F. J. Zerilh, Phys. Rev. D 2 (1970), 2141. 

S. A. Teukolsky, Astrophys. J. 185 (1973), 635. 

S. Chandrasekhar, Mathematical Theory of Black Holes (Oxford University Press, 1983). 
Y. Mino, M. Sasaki, M. Shibata, H. Tagoshi and T. Tanaka, Prog. Theor. Phys. SuppL 
No.128 (1997), 1; gr-qc/9712057 

M. Sasaki and H. Tagoshi, Living Rev. ReL (2003), to appear.; |^qc/0306120l 
B. S. DeWitt and R. W. Brehme, Annals Phys. 9 (I960), 220. 



W. Hikida et al. 



15) Y. Mino, M. Sasaki and T. Tanaka, Phys. Rev. D 55 (1997), 3457; |gr-qc/9606018| 
Prog. Theor. Phys. Suppl. 128 (1997), 373; gr-qc/9712056 
T. C. Quinn and R. M. Wald, Phys. Rev. D 56 (1997), 3381; gr-qc/9610053 

S. Detweiler and B. F. Whiting, Phys. Rev. D 67 (2003), 024025; gr-qc/0202086 

S. Mano, H. Suzuki and E. Takasugi, Prog. Theor. Phys. 95 (1996), 1079; |gr-qc/9603020| 

Prog. Theor. Phys. 96 (1996), 549; gr-qc/9605057 ^ 

S. Mano and E. Takasugi, Prog. Theor. Phys. 97 (1997), 213; |gr-qc/9611014| 
P. L. Chrzanowski, Phys. Rev. D 11 (1975), 2042. 
A. Ori, Phys. Rev. D 67 (2003), 124010; gr-qc/0207045 

Y. Mino, H. Nakano and M. Sasaki, Prog. Theor. Phys. 108 (2003), 1039; 'gr- qc/0111074| 
L. Barack and A. Ori, Phys. Rev. D 61 (2000), 061502(R); gr-qc/9912010 
H. Nakano, Y. Mino and M. Sasaki, Prog. Theor. Phys. 106 (2001), 339; gr-qc/0104012 
L. Barack, Y. Mino, H. Nakano, A. Ori and M. Sasaki, Phys. Rev. Lett. 88 (2002), 091101; 
gr-qc/0111001 

L. Barack and A. Ori, Phys. Rev. D 66 (2002), 084022; 'gr-qc/0204093l 
L. Barack and A. Ori, Phys. Rev. D 67 (2003), 024029; gr-qc/0209072 
L. Barack and A. Ori, Phys. Rev. Lett. 90 (2003), 111101; gr-qc/0212103'' 
S. Detweiler, E. Messaritaki and B. F. Whiting, Phys. Rev. D 67 (2003), 104016; 
gr-qc/0205079 

L. M. Burko, Phys. Rev. Lett. 84 (2000), 4529; gr-qc/0003074' 

N. Sago, H. Nakano and M. Sasaki, Phys. Rev. D 67 (2003), 104017; gr-qc/0208060j 
H. Nakano, N. Sago and M. Sasaki,; gr-qc/0308027 
Y. Mino, Phys. Rev. D 67 (2003), 084027; gr-qc/0302075 

P. Jaranowski and G. Schafer, Phys. Rev. D 57 (1998), 7274 [Errata 63 (2001), 029902]; 
gr-qc/9712075 

P. Grandclement, E. Gourgoulhon and S. Bonazzola, Phys. Rev. D 65 (2002), 044021; 
gr-qc/0106016 

T. Damour, B. R. Iyer and B. S. Sathyaprakash, Phys. Rev. D 57 (1998), 885; 
.gr-qc/9708034. 



